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ABSTRACT 

The properties of the Galactic H I shell GSH 138—01—94 are derived 
from data of the Canadian Galactic Plane Survey. The basic parameters of 
GSH 138—01—94 were determined by fitting the expansion of a thin shell to 
the expansion velocity field on the sky. The kinematic distance is 16.6 kpc for 
"i^LSR = —94.2 ± 0.5 km s^^ The radius is 180 ± 10 pc, the expansion veloc- 
ity fexp = 11-8 ± 0.9 km s~^, and the mass 2 x 10^ Mq. No radio continuum 
counterpart of the shell was detected at 21 cm or at 74 cm. Absorption of a 
background continuum source constrains the spin temperature of H I in the shell 
to Ts = 230?^| K. The expansion age of GSH 138-01-94 is 4.3 Myr. These ob- 
servables are in excellent agreement with predictions from hydrodynamic models 
for a supernova remnant in a low-density low-metallicity environment such as 
the outer Galaxy. GSH 138—01—94 is then the largest and the oldest supernova 
remnant known. It provides direct evidence for the release of mechanical energy 
in the interstellar medium by stars in the outer galaxy. It is argued that such old 
supernova remnants be found in low-density, low-metallicity environments such 
as the outer Galaxy, dwarf galaxies and low surface brightness galaxies. 

Subject headings: ISM: bubbles, kinematics and dynamics, supernova remnants. 
Galaxy: kinematics and dynamics 
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1. Introduction 

The rate and the location of star formation is an important factor in the energy budget 
of the interstellar medium and, in the long term, the evolution of galaxies. The gaseous disk 
in most galaxies is much more extended than the stellar disk, leaving a significant fraction of 
the mass of the interstellar medium at large radii outside the bright optical disk. Despite the 
significant gaseous mass at large radii, little or no star formation occurs outside the stellar 
disk. 

The most widely accessible tracer of kinetic energy in the interstellar medium inside 
and outside the stellar disk of galaxies is the H I velocity dispersion, am- Remarkably, am 
remains constant, approximately 10 km s~^, with distance from the centre, despite a drop of 
more than two orders of magnitude in the surface density of stars (Van Zee & Bryant 1999). 

Energy sources that have been proposed to sustain the ~ 10 km s~^ H I velocity dispersion 
in the outskirts of galaxies include previously undetected low-level star formation (Ferguson 
et al. 1998a), magneto- hydrodynamic instabilities in the disk (Sellwood & Balbus 1999), and 
infalling gas clouds (Hunter et al. 1998, 1999). Which of these, if any, is the dominant source 
of energy is still an open question. 

In this paper we present evidence for supernova activity in the extreme outer region of 
our Galaxy. Evidence for star formation at large distances from the Galactic centre exists 
from detection of molecular gas, distant O and B stars, and their associated H II regions. In 
this paper, we present evidence for an expanding H I shell in the same area of the Galaxy as 
surveyed by Digel et al. (1994). The structure is visible in the Maryland-Greenbank survey 
map shown by Digel et al. (1994), but the southern half is so faint that it was not previously 
recognized as a shell. 

2. Observations 

2.1. The Canadian Galactic Plane Survey 

The data presented in this paper were obtained as part of the Canadian Galactic Plane 
Survey. H I 21-cm line observations were carried out with the Dominion Radio Astrophysical 
Observatory (DRAG) interferometer, supplemented with zero-spacing information from the 
DRAG 26m telescope. These data represent all structure down to the 1' x 1'/ cos 5 resolution 

Fig. 1. — H I channel maps from the DRAG 21-cm line data. Every other channel between 
—78 km s~^ and —110 km s~^ is shown. 
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limit of the DRAO interferometer at 21 cm wavelength and dechnation 5. Additional obser- 
vations of 21-cm continuum and 74-cm continuum from DRAO, and CO spectral line data 
from the Five College Radio Astronomy Observatory (FCRAO) were used. A description of 
these data sets can be found elsewhere (Higgs et al. 1999; Taylor et al. 1999; Heyer et al. 
1998; Carpenter 1997). 

3. Results 

3.1. General morphology 

Figure 1 shows the shell that we refer to as GSH 138—01—94 , as a circular structure 
with radius 37f3, best visible in the channel maps around velocity —96 km s~^. The blue- 
shifted side of the shell represents the most extreme velocity of all H I in this direction. 
On the redshifted side of the shell, the confusion with the background is stronger. The 
channel maps between —88.0 km s~^ and —102.9 km display a distinct empty region 
surrounded by areas of higher intensity. The redshifted half of the shell is clearly fainter 
than the blueshifted half. However, the redshifted half shows several filamentary structures 
similar to blueshifted filaments that can be safely attributed to the shell. Examples of 
such features are found at {l,b,v) = (138?3, -1?8, -89.7) and (138?0,-0?7,-88.0). At an 
even more redshifted velocity the empty region inside the ring fills up to become a local 
maximum, most clearly visible in the channel maps of —81.4 km s~^ and —83.1 km s~^, that 
disappears as the velocity decreases to —78 km s~^. On the blueshifted side, the ring shape 
starts fiUing up in the channel maps at —106.2 km s~^. Therefore, the change in morphology 
with velocity is symmetric with respect to the channel map at —94.6 km s~^, although the 
redshifted side is somewhat fainter than the blueshifted side. This behavior is characteristic 
of an expanding shell, with both the near (blueshifted) and far (redshifted) sides visible. At 
the lowest intensity levels visible in Figure 1, the shell is nearly continuous, but the side that 
faces the Galactic equator is a factor 2 brighter than the opposite side, on average. 

The shell is further illustrated in the latitude- velocity diagram in Figure 2. The red- 
shifted emission that we associate with GSH 138—01—94 traces the doppler ellipse very 

Fig. 2. — Latitude-i'LSR diagram, constructed by adding three latitude-velocity slices at 
longitude / = 137?865 and 3' offset in longitude on either side. At this longitude, which 
is slightly offset with respect to the centre of the shell, the slices intersect gas at the most 
extreme velocity. The ellipse represents the thin shell model fitted in section 3.1, taking into 
account that the shell is not intersected through the centre in this Figure. 
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closely. Furthermore, this emission, with velocity —83 km s""^ is restricted to the latitude 
6 ~ — 1?25. The emission is therefore not associated with a more widely distributed compo- 
nent, like most of the emission in Figure 2. 

The centre, radius, expansion velocity and central velocity of the shell were determined 
by fitting the velocity field of a thin expanding shell to the data at a number of positions. The 
line of sight velocity ^^los for a thin expanding shell with centre {xcUc), expansion velocity 
Vq, radius Rs, and central velocity Vc, at any point {x,y), with {x — x^Y + {y — VcY < 
is given by 



The five parameters defining this function can be determined by fitting a large number of 
measured points (xj, t'j). In particular, the expansion velocity of the shell is determined 
by on and off-axis data rather than only by emission at the extreme velocities. Similarly, 
the radius is determined from points over many channels rather than from a single channel 
map. 

The data points {xi,yi,Vi) may in principle be distributed anywhere within the radius 
of the shell, although in practice, there are some restrictions. Measurements well within 
the outer boundary are required to allow a determination of the expansion velocity. Mea- 
surements near the outer boundary of the shell should be avoided, because the component 
of the expansion velocity along the line of sight vanishes there. Furthermore a degeneracy 
between the radius and expansion velocity of the shell exists if all measured points are at 
the same distance from the centre. Data from both sides of the shell are not required, but 
more accurate results can be expected if data from both the approaching and receding side 
are included. 

Line of sight velocities at 94 positions throughout the area of GSH 138—01—94 were 
determined by fitting the local line profile with 3 to 5 Gaussians and a second degree poly- 
nomial baseline. We are only concerned with obtaining a good fit to the centroid of the 
one component in the line profile associated with the shell. The remaining components were 
added to achieve a satisfactory fit to the H I line profile. Only the position and velocity of the 
Gaussian component identified with the shell were used. The distribution of the measured 
points on the sky is shown in Figure 3. More data exist for the approaching side of the shell 
because this side does not suffer from confusion with other emission. The results of the fits 
(see Table 1) did not change by a large amount when the fit was repeated for five randomly 
selected subsets containing only one third of the data. The spread in the results from the 
fits to subsets of the data was used to determine the errors listed in Table 1. The results 
are consistent within the errors with estimates of the geometric centre and radius defined 




(1) 
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Fig. 3. — Left: Locations of data points used for the fit on the sky. Symbols: triangles: 
\v\os — Vc\ > 10 km s^"*^, squares: 10 km > \vios — Vc\ > 5 km s^"*^, circles: 5 km s^^ > 
|i^ios — '*^c| > km s~^. Filled symbols represent the receding (far) side of the shell, open 
symbols the approaching (near) side. Right: Line of sight velocity of shell material as a 
function of projected distance from the centre. The curved hues are the best fit solution of 
expression (1) (solid line) and the solutions with la perturbations of radius and expansion 
velocity (dashed and dot-dashed curves). 



Table 1: Properties of GSH 138-01-91 . 



Quantity 


Value 


Reference 


Centre (Galactic) 


I = 138°6' ± 3' 6 = -1°15:5 ± 3' 


Section 3.1 


Centre (J2000) 


a = 2'^50™35!0 S = +58°00'15" 




Central velocity 


-94.2 ± 0.5 km s'^ 


Section 3.1 


Expansion velocity 


11.8 ±0.9 km s-i 


Section 3.1 


Heliocentric distance 


16.6 ±2.5 kpc 


Section 3.2 


Galactocentric distance 


24 kpc 


Section 3.2 


Radius 


37:3 ± 2.2 


Section 3.1 




180 ± 10 pc 




Mass 


2 X 10^ M0 


Section 3.6 


Expansion kinetic energy 


3 X 10^3 J 


Section 3.6 


Expansion momentum 


2 X 10^ Mo km s"^ 


Section 3.6 


Expansion age 


4.3 X 10^ year 


Section 3.6 
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by the average and half the difference of the coordinates of opposite sides of the shell in the 
channel maps. 

Adopting the centre derived from the fit, we plot 

^los ^ 

function of projected distance 

from the centre in order to display the data and the fit in a 2-dimensional graph. This graph 
is shown in Figure 3. The data points trace the Doppler ellipse outhned in Figure 2. The 
scatter in the data exceeds the formal error in the Gaussian fits, and likely reflects departures 
from the idealized expansion law assumed in equation (1). The solid line represents the best 
fit to all the data, and the interrupted lines indicate the errors quoted for the radius and 
expansion velocity. Note that a smaller expansion velocity requires a larger radius to obtain 
a reasonable fit. The most extreme negative velocities are not very well represented by the 
flt. These points can be traced to the vertical fllament visible in the channel maps around 
-107.8 km s-\ 



3.2. Distance 

The central velocity of GSH 138—01—94 is close to the extreme velocity of Galactic H I 
. This implies a large galactocentric radius for any choice of Galactic rotation curve that 
is reasonably flat. As many derived properties depend strongly on the adopted distance, a 
careful review of the uncertainties involved is essential. 

The line-of-sight velocity of a cloud on a circular orbit with radius R and velocity V{R) 
around the Galactic centre assuming circular orbits, is 

^ios= (^n^)-l^(i?o))sinZ 

Here, Rq is the distance of the sun to the Galactic centre, and / is the Galactic longitude 
of the cloud. This relation can be solved for the Galactocentric distance R, and hence the 
heliocentric kinematic distance d. For GSH 138—01—94 , = —94.2 km s^^, and / = 138°6', 
so R = 23.6 kpc, and d = 16.6 kpc, assuming V{R) = 220 km s~^ and Rq = 8.5 kpc. 

If GSH 138—01—94 is not at such a large distance, then it must have a high anomalous 
velocity with respect to the surrounding medium. We rule this out for two reasons. First, 
the appearance of GSH 138—01—94 on the sky is nearly perfectly circular. This effectively 
rules out a velocity relative to its immediate surroundings which is a large fraction of the 
expansion velocity (11.8 km s"^). Second, the mass of GSH 138—01—94 is hkely dominated 
by swept-up interstellar gas. Suppose the source of the shell has a velocity of 100 km s~^ 
relative to the ambient medium, and the mass of the ejecta amounts to 10% of the mass of 
the shell. Conservation of the momentum of the ejecta due to the velocity of the source. 
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implies that the pecuhar velocity of the shell would be only 10 km s ^ in this example. 

Of course, this argument does not exclude that the velocity of the surrounding medium 
is also different from standard Galactic rotation. We now consider locahzed departures from 
Galactic rotation along the line of sight, which could create the false impression of a large 
distance, and departures from standard Galactic rotation in the outer Galaxy. 

Large-scale streaming motions in galaxies are associated with bars or spiral arms. 
Streaming motions associated with spiral arms are the only relevant factor outside the so- 
lar circle. The spiral shock model of Roberts (1972) for the longitude range I — 130° to 
I — 140° allows a minimum line of sight velocity of approximately —55 km s~^, which is some 
20 km s~^ blucshiftcd from standard Galactic rotation. The velocity of GSH 138—01—94 is 
not allowed by the density wave model for the Perseus arm by a margin of 40 km s^^. Also, 
absorption of radio continuum emission towards the W3/W4 H II regions is not observed 
beyond t'LSR = —50 km (Normandeau 1999). We therefore conclude that the velocity of 
GSH 138—01—94 implies that it is located well beyond the Perseus arm. 

The validity of standard Galactic rotation at large radii is uncertain. The rotation 
velocity of most spiral galaxies remains nearly constant to the outermost point that can be 
measured, and the same is true for the Galaxy to a distance 2i?o from the centre (Merrifield 
1992). By comparison, the kinematic distance of GSH 138—01—94 imphes it is 2.8i?o from 
the Galactic centre. The issue concerning the distance to GSH 138—01—94 is not whether the 
rotation curve of the Galaxy declines at large distances, but rather whether the assumption 
that the gas follows circular orbits is vahd. The magnitude of this effect may be estimated 
from the / — > — / asymmetry of the extreme velocity of H I in the Galaxy (Kuijken Sz Tremaine 
1994). The magnitude of this asymmetry is at most 20 km s^^, and the sign implies a more 
negative velocity in the second quadrant. We estimate the effect of non-circular orbits by 
adding a radial expansion term of —10 km s~^ to the velocity V{R) and recalculating the 
kinematic distance. The result is a distance which is 3 kpc or 20% smaller than the distance 
derived assuming circular orbits. 

In this paper we therefore adopt a distance of 16.6 kpc for GSH 138—01—94 , keeping 
in mind that a systematic effect of the order of 20% may exist due to non-circular orbits in 
the outer Galaxy. 

3.3. Environment 

The large distance of GSH 138—01—94 places it in the far outer Galaxy. The Galactic 
latitude of its centre corresponds with a distance 370 pc south of the Galactic plane. The 
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actual centre of the H I layer cannot be determined with a high accuracy for a particular 
longitude. From Figure 60b in Burton (1991) it is estimated that the centre of the H I layer is 
between 0.5 and 1 kpc toward positive Galactic latitude in the vicinity of GSH 138—01—94 . 
The half-width thickness of the Galactic H I layer at 20 kpc from the centre is approximately 
1.5 kpc (Merrifield 1992). Therefore, GSH 138—01—94 is not exceptionally far from the 
local midplane, and its diameter is about a quarter of the thickness of the disk at half 
maximum. The side of GSH 138—01—94 in the direction of the Galactic equator is a factor 
two brighter than the side facing the opposite direction, which is consistent with a density 
gradient perpendicular to the plane. 

Adjacent to GSH 138—01—94 at lower longitude an H I complex is observed that may 
have some connection with GSH 138—01—94 . The H I complex displays a pecuhar velocity 
gradient in the sense that the velocity is more negative at a larger distance from the Galactic 
equator. For a disk in which the rotation speed does not change with distance from the 
midplane, the velocity is expected to become less negative by a factor cos (6), which is 
negligible over the latitude range considered. Although this H I complex may be physically 
related to GSH 138—01—94 , it may also be a coincidence because a small velocity interval 
corresponds with a large distance along the line of sight in this part of the Galaxy. 

The gas density in the vicinity of GSH 138—01—94 was derived according to the relation 
nn = Nni/Ad), with Nhi the H I column density obtained by integrating the brightness tem- 
perature over a narrow velocity interval Av centered on the central velocity of the shell. Ad 
is the length of the line of sight corresponding to Av and Galactic rotation. The brightness 
temperature averaged over an area that covers the latitude range of GSH 138—01—94 on 
both sides in longitude is — 13 K, and does not vary significantly with velocity for 
Av < 10 km s~^. The resulting hydrogen density is riu — 0.02 cm~^, which is fairly inde- 
pendent of the particular choice of Av. This density is within a factor 2 of the azimuthally 
averaged H I density riu — 0.04 at R = 2.35i?o in Figure 8 of Merrifield (1992). A conse- 
quence of the steep slope of the velocity-distance relation for the longitude and velocity of 
GSH 138—01—94 is that the derived density is the average across a significant line of sight: 
Ad = 4.7 kpc for Av = 10 km s~^. This will be discussed further in Section 3.6. 



3.4. Substructure 

At the resolution of 1 arcminute a great deal of substructure is observed. At the most 
extreme velocity (i'lsr — —110 km s""*^) a prominent filament is seen in the channel maps. 
Such condensations and filaments may form as a result of the Rayleigh- Taylor instability 
if a dense shell is expanding due to the pressure of the tenuous hot interior, and thermal 
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instabilities because the cooling increases with density. A number of compact spots can be 
seen in Figure 1, e.g. at {l,b,v) = (138?125, -0?695, -106.2). 

At (Z, 6, v) = (137?79, — 0?81, —99.6) and adjacent channels, a conspicuously empty, 7'.2 x 
4f2 eUiptical area is visible in the rim of the shell. Over the velocity —90 km s~^ > vlsr > 
—110 km s~^, the brightness temperature remains below 15 K in this area. This is consistent 
with the intensity of the surroundings of the shell, and about half the intensity of the bright 
edge surrounding this area. The isophotes of the hole are elliptical, with major axis in 
position angle 36°, and axial ratio 0.59. The projected distance of the hole from the centre 
of GSH 138—01—94 is 0.895 times the radius of the shell in position angle —33°. A circular 
hole in the shell seen in projection at the same location, would appear as an elliptical hole 
with its major axis in position angle 33°, and axial ratio 0.45. Taking into account the 
uncertainties in the positions, the observed shape of the hole is consistent with a nearly 
circular hole in GSH 138—01—94 . At the distance of 16.6 kpc, the diameter of this hole is 
35 pc. We refer to this hole as Cavity 1. 

A smaller, but very distinct empty region is observed at (Z, 6, v) = (137?54, — 1?05, —101.2), 
and the adjacent channel maps. The central brightness temperature remains consistent with 
the background, which is a factor 3 below its immediate surroundings. We refer to this 
feature as Cavity 2. 

3.5. Spin temperature 

A weak but distinct absorption feature is observed against a, T^, — 82.8 K background 
continuum source in the direction (Z,6) = (137?93, — 1?27). Figure 6 A shows the absorp- 
tion profile in this direction. We determine the peak brightness temperature by fitting a 
Gaussian and a constant offset to the data in a specific velocity interval. The sum of the 
Gaussian and the offset is a relatively accurate estimate of the brightness temperature at 
Vlsr = —107 km s^^, despite the modest signal to noise ratio of typically 7, because it uses 
information from more than one channel. This was done for eight positions, ±3' offset in 
longitude or latitude, from the position of the continuum source. The fitted offset in intensity 
contributes typically less than 20% to the total intensity. The mean brightness temperature 
at the eight off-positions is Tb^oflf — 18-6 ± 1.3 K. The error indicates the rms variations, 
divided by -\/8, and is comparable to the formal error in the fit of a single spectrum. At the 
position of the continuum source the brightness temperature is 7b,on = 11-9 ± 1.5 K, which 
is less than any of the eight surrounding positions. These numbers imply a spin temperature 
between 173 K and 367 K, for la variation of Tb^on and T\^^oS, with a most likely value of 
Tg = 230 K. The optical depth is 0.084. Although velocity crowding can place unrelated 
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gas at the same velocity, the fact that we see a narrow absorption feature at this velocity 
argues against this possibility. We therefore conclude that the spin temperature derived here 
applies to H I in GSH 138-01-94 . 

An upper limit to the kinetic temperature can be derived from the width of the line 
profiles. The average H I velocity dispersion is 1.9 km s~^ for the eight off-positions, with 
scatter consistent with the formal error of 0.4 km s~^ from the fit. Fits at other locations 
in the shell give a similar result. An upper limit Tkin < 430 K to the kinetic temperature is 
found by assuming that the observed linewidth is completely due to thermal broadening. This 
upper limit is consistent with the low spin temperature derived above, and emphasizes that 
the temperature in the shell is not far from the temperature of warm H I in the interstellar 
medium. 



3.6. Mciss, energy, momentum 

The relatively smooth H I background visible in Figure 1 was subtracted by applying 
a median filtering procedure outlined in Section 3.7 to the H I channel maps, with one 
difference. Instead of applying a cutoff in intensity to exclude bright point sources, the 
brightest parts of the shell were masked manually in each channel map. Inspection of the 
fine profiles next to the brighter parts of GSH 138—01—94 showed that the zero level in 
the residual data cube is well defined. Some uncertainty is unavoidable because emission of 
GSH 138-01-94 is blended with the H I cloud on the right hand side of GSH 138-01-94 in 
the velocity interval —94.6 km s~^ to —101.2 km s~^, around longitude 137?5. 

Summation of the background subtracted data cube between —82.3 km s~^ and —115.2 km s 
results in a total H I mass Mm — 1-8 x 10^ within a circular area defined by the out- 
line of GSH 138-01-94 . The area bounded by the limits / < 137?75, b > -1?45, and 
the outer radius of GSH 138—01—94 is affected by confusion with the adjacent cloud (see 
Figure 1). The H I mass in this area is 4.2 x 10"^ M©, or 23% of the total. In order to 
correct for the presence of the cloud, we assume that the H I mass in the area affected by 
confusion is equally distributed between the shell and the cloud. With this correction, we 
find Mhi = 1.6 X 10^ Mq. The atomic hydrogen mass is an order of magnitude larger than 
the molecular mass discussed in Section 3.8. Taking into account primordial Helium, we 
find the total mass of GSH 138—01—94 is Mtot = 1.3Mhi. The kinetic energy associated 
with the expansion of the shell then follows from E'kin = ^MtotV^xp the momentum from 
Pexp = MtotVexp- Numerical values for these quantities are listed in Table 1. 

The hydrogen mass in the volume of the shell assuming the density riu — 0.02 cm~^ 
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(Section 3.3) is 1.2 x 10^ M©, an order of magnitude smaller than the H I mass of the 
shell. This difference cannot be the result of an error in the distance. The swept-up mass 
is proportional to the volume, therefore the third power of distance. The HI mass derived 
from the linefiux is proportional to the square of the distance. Therefore, a ten times larger 
distance would have to be assumed to obtain consistent results. The discrepancy can be 
resolved if the environment of GSH 138—01—94 has a higher density than the Galactic 
average at the Galactocentric radius of GSH 138—01—94 . The density found in Section 3.3 
is an average over a long line of sight. If the swept-up mass and the mass of the shell are 
equal, the density is = 0.27 cm""^, seven times the average midplane density found by 
Merrifield (1992). This estimate of the ambient density is preferred over the density derived 
in Section 3.3, but the uncertainty in the density is at least a factor 2. 

The expansion age of GSH 138—01—94 is Texp = aR/vexp, assuming the expansion law 
R t"'. We assume a — ^. This is apphcable to a shell expanding from the pressure of a 
hot interior, such as an old supernova remnant (McKee & Ostriker 1977). The expansion 
age thus becomes Tcxp = 4.3 Myr. The precise value of a is of less importance given the 
uncertainties, but we note that we somewhat underestimate the age in the likely event that 
the expansion law was steeper in the past. 

3.7. Radio continuum 

Figure 4 shows the DRAO 74-cm continuum map of the area of GSH 138—01—94 . The 
gradient in the background intensity in Figure 4 is due to the Galactic background emission 
that decreases with distance from the Galactic equator. There is no indication of enhanced 
emission, on top of the Galactic background, associated with GSH 138—01—94 . In order to 
derive an upper limit to the continuum flux of GSH 138—01—94 , we subtracted the smooth 
Galactic background emission in the following way. 

First, point sources brighter than approximately Tb = 20 K above the diffuse background 
were excluded from the analysis by clipping on intensity. The clipped image was subjected 

Fig. 4. — Gray scale representation of the 408 MHz continuum image. Gray scales are linear 
and range from 50 K (white) to 65 K (black). The location of the H I shell is indicated 
by the white contour, of the column density 2 x 10^° cm~^, obtained by integration over 
the velocity range —98 km s~^ to —107 km s~^. A + marks the expansion centre of 
GSH 138—01—94 determined in Section 3.1. The square indicates the field of view of Figure 5. 
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to a spatial median filter with dimensions 6' x 90' in latitude and longitude, in two steps. 
In the first step, the median filtered image was subtracted from the original image in order 
to obtain a first order residual image of structure on scales much smaller than the size of 
the median filter. This first-order residual image was clipped at the Tb = 4 K level, setting 
all pixels below this level to 0. Only individual point sources exceeded the 4 K level. The 
clipped first-order residual image was then subtracted from the original image. 

In the second step, the median filter was applied to the point source subtracted image. 
This second pass resulted in our final estimate of the smooth background. The residual 
image after subtraction of the final smooth background image, has a well-defined zero level. 
The elongated shape of the median filter allows effective subtraction of structure that is 
extended in Galactic longitude, while retaining elongated structure oriented at a large angle 
with respect to the Galactic plane. Negative residuals are less than 7% (typically 2%) of 
the smooth background emission in the vicinity of GSH 138—01—94 . The background 
subtraction procedure was tested in another area of the CGPS, where continuum shells of 
the same angular size as GSH 138—01—94 were detected. It was confirmed that a shell 
structure of the size of GSH 138—01—94 is effectively separated from the background by this 
procedure. This method, then, effectively allows us to identify a shell in continuum, if it 
exists. No shell was observed above the noise. 

The sensitivity of the DRAG 408 MHz map is limited by confusion of faint background 
sources. Excluding individually visible point sources, the r.m.s. fluctuations in the residual 
image arc ATb = 0.86 K. Due to the crowding of the many resolved point sources, it is not 
feasible to integrate the flux over the area of GSH 138—01—94 . Instead we adopt the r.m.s. 
fluctuations of the background as a safe upper hmit to the surface brightness. An object 
with the size of GSH 138—01—94 and brightness temperature Tt, = 0.86 K (Icr) would have 
been recognized in the residual map without difficulty. The number of independent beams 
over the area of GSH 138—01—94 should not be used to quantify the statistical significance 
of this upper limit, because it does not take into account the many points sources in the 

Fig. 5. — Detail of GSH 138—01—94 displaying the area around Gavity 1 indicated by the 
box in Figure 4. Gray scales represent H I column density, integrated over the velocity 

range —98 km s~^ to —107 km s~^. The range of gray scales is 0.6 x 10^'^ cm^^ (white) to 
8.0 X 10^'^ cm~^ (black). White contours give the CO brightness temperature integrated over 
the same velocity range. Contour levels for the CO map are 2 (Sex), 3, 4, 5 K km s^^. The 
21-cm continuum intensity is shown as black contours at T], = 0.14 (2(t), 0.28, 0.42, ... K 
above the mean background level of 5.13 K. The -|- indicates the position of the Bl star 
studied by De Geus et al. (1993). The peaks in the Ha emission are indicated with x. 
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field. 

The 408 MHz surface brightness is then bounded by S408 < 4 x 10~^^ W m"^ Hz~^ sr"^ 
This upper hmit is significant, because it is an order of magnitude below the surface bright- 
ness of known old supernova remnants (Landecker et al. 1990). An upper hmit to the 408 
MHz flux density can be derived from the upper limit to the surface brightness. There is 
some uncertainty in the choice of the area to which the surface brightness limit applies. A 
likely possibility is that a radio continuum counterpart of GSH 138—01—94 would have a 
shell morphology similar to the HI. In this case, the upper limit to the surface brightness 
is really a limit to the surface brightness of the bright outer rim of such a continuum shell. 
Assuming a ring with outer radius 37^3, inner radius 30', and a uniform surface brightness 
equal to the upper limit translates in an upper limit of 0.57 Jy to the flux density at 408 
MHz. 

Figure 5 shows the 1420 MHz continuum image of the area indicated by the square in 
Figure 4. An extended continuum source is seen in projection on the edge of Cavity 1. When 
rotated by 20 degrees counter clockwise. Figure 5 can be compared directly with Figure 3 
in De Geus et al. (1993). Radio continuum emission is observed at the location of the two 
maxima in Ha intensity in the map of De Geus et al. (1993), as indicated in Figure 5. In fact, 
the H II region is extended and covers the same area on the sky as the faint resolved 21-cm 
continuum emission. Based on the correlation with Ha and the absence of a counterpart 
at 74 cm, the continuum emission may well be thermal emission from the H II region. The 
velocity of the H II region as determined by De Geus et al. (1993), is —101 km s~^. The H I 
that can be associated with this H II region and the adjacent molecular cloud (Section 3.8) 
is part of the edge of Gavity 1, which is a hole in GSH 138—01—94 (Section 3.4). 

3.8. MoleculEir gas 

The FCRAO CO data show a small molecular cloud (number 2 in Digel et al. 1994) 
at {l,b,v) = (137?76, —0?97, —103.6), with peak brightness temperature Tb = 2.9 K, and 
velocity dispersion aco = 0.94 km s^^. De Geus et al. (1993) found an H II region with the 
same velocity as the molecular cloud and a candidate ionizing B star with a photometric 
distance between 17.5 and 32 kpc. The relative locations of the molecular cloud, the H II 
region and the ionizing star are indicated in Figure 5. 

Figure 6 B shows the CO and H I spectra at this position. H I emission at vlsr ~ 
— 101.9 km s~^ is associated with GSH 138—01—94 . The H I and CO proflles can be com- 
pared directly thanks to the arcminute resolution and zero spacing information in the DRAO 
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Fig. 6. — A: Absorption of Galactic H I against the continuum source at (/, b) = 
(137?93, — 1?27). The smooth emission spectrum is the average of eight positions 3' off- 
set in I and b. The negative curve is the difference between the spectrum at the position of 
the continuum source and the average spectrum of the surroundings. Note the absorption 
feature at vi^sr — —108 km s~^, which is associated with GSH 138—01—94 . B: CO emission 
hne (histogram) and H I emission in the direction of the CO cloud. 
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H I data. The coincidence of these two features in position and velocity is highly unlikely 

to be a chance alignment. We therefore associate the CO cloud with the approaching side 
of the H I shell. Note that this association of the CO cloud with GSH 138—01—94 means 
its distance is 16.6 kpc (Section 3.2), considerably smaller than the kinematic distance of 
21 kpc adopted by Digel et al. (1994). It is noteworthy that cloud 2 is the only one of 11 
clouds in the sample of Digel et al. (1994) for which the mass derived from the CO flux ex- 
ceeds the virial mass. The smaller distance that we derived here, decreases this discrepancy 
significantly, because the virial mass scales only linearly with distance. 

In this interpretation, this molecular cloud is only marginally more distant than the 
other molecular clouds in the sample of Digel et al. (1994), and only marginally inconsistent 
with the photometric distance of the B star found by De Geus et al. (1993). Digel et al. 
(1994) found the mass of this molecular cloud to be 2.3 x 10^ Mq, assuming a distance of 
16.6 kpc. The bright northern rim of the GSH 138—01—94 was targeted by Digel et al. 
(1994), but no CO emission was found there. Digel et al. (1994) identified the pointsource 
IRAS 02450+5816 with their cloud 2, but there is no indication of a shell structure in the 
IRAS maps. 



4. Discussion 

4.1. Source of the shell 

The kinematic distance to GSH 138—01—94 is the basis for our interpretation of the data 
presented in this paper. Summarizing the arguments given in Section 3.2, we conclude that 
systematic effects between 10 km s^^ and 20 km s^^ may exist, but these do not seriously 
affect our conclusions. We estimate an error in the distance of 20%, similar to the assessment 
of Digel et al. (1994) for their molecular clouds. The association of the H II region found 
by De Geus et al. (1993) and a photometric distance of its probable ionizing star with 
GSH 138—01—94 supports the large distance derived in Section 3.2. The conclusion that 
GSH 138—01—94 is in the far outer Galaxy does not depend on the exact shape of the 
rotation curve in the outer Galaxy. The distance of GSH 138—01—94 is therefore sufficiently 
secure to allow a meaningful comparison with other shells of known origin. 

The size of GSH 138—01—94 is similar to small superbubbles, such as the Orion-Eridanus 
bubble (Brown et al. 1995). However, the kinetic energy of the expanding shell (3 x 10"^^ J) 
is consistent with the lO"*"^ J released in a single supernova explosion (Leitherer et al. 1992). 
Although the expansion energy of GSH 138—01—94 is a factor 3 smaller, this is well within 
the uncertainty in the energy release of a supernova, and which fraction of this energy is in 
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the form of kinetic energy of the ejecta. A recent discussion was given in Thornton et al. 
(1998). The energy released in a stellar wind can be of the same order of magnitude as the 
kinetic energy of the ejecta of a supernova explosion. Wc first explore what star could be 
the source of GSH 138—01—94 assuming the shell is the result of only a stellar wind. 

Normandeau et al. (2000) discuss the possible sources of the H I shell G132.6— 0.7— 25.3 

(radius 33 pc). This shell may be a bubble blown by the stellar wind of the Bl la star 
BD60°447, but the expansion age of this shell is less than the age of the star. Normandeau 
et al. (2000) conclude that G132.6— 0.7— 25.3 can also be an old supernova remnant. Brans- 
ford et al. (1999) found shells with radii up to 80 pc in a survey of Wolf-Rayct stars in M31. 
Although GSH 138—01—94 is more than twice the size of the largest of these shells, it cannot 
be excluded that this is the result of a lower ambient density. The most significant param- 
eters for the interpretation of GSH 138—01—94 as a wind blown bubble are the expansion 
kinetic energy and the age of the shell. The expansion kinetic energy of GSH 138—01—94 is 
30 times larger than that of G132.6— 0.7— 25.3 (Normandeau et al. 2000). Note that Norman- 
deau et al. (2000) found only one star that could have been the origin of G132.6— 0.7— 25.3 
by means of its stellar wind. 

The mass loss rate from a stellar wind increases rapidly with luminosity. The longer 
lifetime of a star with lower mass compensates only partly for the much smaller mass loss 
rate. The mass loss rate also depends on the metallicity Z. Nugis & Lamers (2000) find that 
the mass loss rate scales as Z™, m = 0.5. Vink et al. (2001) predict m = 0.64 to m = 0.69, 
and argue that oxygen is the relevant trace element for Z at low metallicity. At 18 kpc from 
the Galactic center, the abundance of oxygen relative to hydrogen is a factor 10 below the 
solar value (Smartt & RoUestone 1997). The velocity of the stellar wind is also smaller at 
low metaUicity (Nugis & Lamers 2000). We conclude that the energy released in the form 
of a stellar wind decreases by a factor ~3 if the metallicity decreases with a factor 10, for 
stars with the same luminosity. 

The relation between stellar luminosity and the energy released per second in the stellar 
wind was defined from the sample of O and B stars of Lamers et al. (1999). The hfetime of 
these stars was related to luminosity at mid-life using the stellar evolution models for solar 
metallicity of Bressan et al. (1993). From this, a relation was found between luminosity and 
the total energy released in the stellar wind during the lifetime of the star. The expansion 
energy of GSH 138—01—94 requires a star with luminosity log(L/L0) = 5.7 at solar metallic- 
ity. Assuming the metallicity is 0.1 solar, we find log(L/LQ) = 5.9. This luminosity implies a 
star with main sequence mass 60 Mq, lifetime 4.1 Myr (Bressan et al. 1993), spectral type on 
the main sequence approximately 05. The same result can be obtained from Abbot (1982). 
It was implicitly assumed that all of the kinetic energy of the stellar wind is available for the 
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expansion of the shell. It is more likely that some fraction is converted into thermal energy 
as the stellar wind and the swept-up interstellar medium pass through a shock. In this case 
a stronger stellar wind (i.e. a more massive star with a shorter lifetime) would be required. 
The age of GSH 138—01—94 is 4.3 Myr assuming an expansion law R ~ t2/7 (Section 3.1), 
similar to the lifetime of a 60 star. It is therefore difficult to provide the expansion 
energy of GSH 138—01—94 in the form of a stellar wind by a star with a lifetime as long as 
the age of GSH 138-01-94 . 

Are the present data consistent with the existence of such a massive star? If the star 
still exists, thermal emission of an H II region associated with the star should be visible in 
the 1420 MHz continuum map, which shows only the small H II region found by De Geus 
et al. (1993) (Figure 5). The effective temperature is lower at a later stage of evolution of 
the star, but only for a brief period of time. In the models of Bressan et al. (1993) the 
effective temperature of a star with initial mass 60 M© is below 3 x 10^ K for only 4000 years 
(0.1% of the lifetime of the star). Assuming an 05V star and the extinction A^/ = 2.7 mag 
found by De Geus et al. (1993), we estimate the star should be magnitude V = 13, observed 
color U — B = —0.5. With this brightness and color, the star would be an obvious target 
for the search for blue stars by Lanning (1977). There is no entry in his list near the center 
of GSH 138-01-94 , although the objects 7 (me ~ 20) and 54 (me ~ 12.7) are within the 
field covered by the channel maps in Figure 1 (but outside the perimeter of the shell) . 

We conchidc that the existence of a star with the mass required by the interpretation of 
GSH 138—01—94 as a stellar wind bubble is unlikely. The radio continuum maps at 408 MHz 
and 1420 MHz also show no evidence for a supernova remnant inside GSH 138—01—94 . If the 
shell is a supernova remnant, the problems raised by the interpretation of GSH 138—01—94 as 
a stellar wind bubble are resolved. In particular, a less massive star is required than in the 
stellar wind case. 

GSH 138—01—94 is larger and older than known supernova remnants. It is illustrative 
to compare GSH 138—01—94 with supernova remnants in the solar neighbourhood. The 
Cygnus Loop has a radius of 19 pc and age 1.4 x 10^ yr (Levenson et al. 1998). The 
Cygnus Loop is therefore of the same size as Cavity 1 (Figure 5). GSH 138—01—94 is larger 
and older than known supernova remnants. Supernova remnants that arc considered to 
be very old include G55.0+0.3 with radius 70 pc, age 1 Myr (Matthews 1998), G65.2+0.6 
with unknown distance but has a low radio surface brightness, (Landecker et al. 1990), 
G106.3+2.7 with radius 70 pc, age 1.3 Myr (Pineauh & Joncas 2000), and G166.2+2.5 
with radius 87 pc, age (Routledge et al. 1986). Among the older remnants would also be 
radio loop I, with radius 60 pc and expansion velocity 25 km s~^ (Heiles 1984) yielding 
an age of 0.7 Myr, assuming this is a single supernova remnant. It is usually difficult to 
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Fig. 7. — Histogram of the number of known supernova remnants as a function of surface 
brightness at 408 MHz. The arrow indicates the upper hmit for GSH 138—01—94 derived in 
Section 3.7. 
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estimate the distance of large shells that can be old supernova remnants, and to show that 
such a shell is the result of a single supernova. Distance ambiguities can be eliminated 
when looking at nearby galaxies. The largest supernova remnant in the optical spectral line 
survey of M31 by Blair et al. (1981) is BA 490, with radius 90 pc. Interestingly, BA 490 is 
located 19 kpc from the centre of M31. Gordon et al. (1999) studied a radio selected sample 
of supernova remnants in M33, eliminating distance ambiguities. The largest supernova 
remnants in their sample have radii of 50 pc, but these authors acknowledge a selection 
effect against the oldest supernova remnants. A direct comparison of GSH 138—01—94 with 
known supernova remnants is therefore inconclusive because cataloged supernova remnants 
are significantly younger, smaller and expand faster. This is likely a selection effect because 
supernova remnants are usually identified in radio or X-ray surveys. As the remnant ages, 
the radio surface brightness declines and the temperature of the interior drops. From this it 
cannot be excluded that an older supernova remnant would be similar to GSH 138—01—94 , 
if it managed to survive long enough. In particular, the fact that no radio continuum shell 
was detected docs not disprove the possibility that GSH 138—01—94 could be a very old 
supernova remnant. 

Figure 7 shows the number of known Galactic supernova remnants as a function of radio 
surface brightness. This histogram is based on the most recent (September 1998) version 
of the catalogue of Green (1984, 1988) available through the CADC. The surface brightness 
at 408 MHz was calculated following Green (1991), adopting a spectral index of 0.5 if no 
spectral index was listed. Six objects out of 220 were excluded because no radio flux was 
listed. The peak of the histogram corresponds with the completeness limit of radio surveys, 
as demonstrated in Figure 3 of Green (1991). 

It is often assumed that an expanding shell dissolves into the ambient medium if its 
expansion velocity becomes similar to the random velocity of clouds, which is of the order 
of 10 km s~^. Glearly this is the case for GSH 138—01—94 , but it is still a well-defined 
shell. The residual velocities in Figure 3 are of the order of 3 km s~^. It is not likely that 
the velocity dispersion in the outer Galaxy is significantly smaller than 10 km s~^. H I 
observations of face-on galaxies show that the H I velocity dispersion is nearly 10 km s~^, 
even far beyond the optical radius (Van Zee & Bryant 1999). The timescale for dissolving 
into the surrounding medium is of the order of t^is = Rs/(^v, where i?s is the radius of the 
shell, and (T„ ^ 10 km is the velocity dispersion of HI clouds. For GSH 138—01—94 we 
find tdis — 18 Myr. McKee & Ostriker (1977) also found that the timescale for the low-density 
cavity created by a supernova remnant will survive approximately 10 times longer than the 
timescale to reach equilibrium between the internal and external pressure. Their tmax uses 
the sound speed in the ambient medium, and is of the same order of magnitude as tdis- The 
long dissolving timescale resolves the apparent contradiction of a well-defined H I shell with 
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expansion velocity of the order of the velocity dispersion of the ambient interstellar medium. 

4.2. GSH 138-01-94 compared with SNR models 

Detailed one-dimensional hydrodynamic simulations of the evolution of supernova rem- 
nants including a realistic, metallicity-dependent, cooling function were made by Thornton 
et al. (1998). A proper treatment of cooling is essential for simulations that proceed beyond 
the adiabatic expansion phase of the remnant. Radiative cooling becomes important first in 
the relatively dense expanding shell. In the late stages evolution of the supernova shell, the 
pressure of the hot interior drives the expansion and compresses the shell further, increasing 
the cooling rate in the shell. This leads to a thin neutral shell in the so-called "pressure 
driven snow plow phase". This qualitative scenario is supported by hydrodynamic simula- 
tions that include radiative cooling, e.g. Cioffi et al. (1988) and Thornton et al. (1998). The 
pressure driven snow plow phase continues until the pressure of the hot bubble drops so far 
that it can no longer drive the expansion of the remnant. Not surprisingly, the timescale 
for the end of the pressure driven snow-plow phase depends on the cooling rate of the gas, 
which itself depends on the density, metallicity and temperature. Thornton et al. (1998) 
included the most detailed cooling function to date and provide relations for the properties 
of the remnant as a function of age, metallicity and density of the ambient medium. 

The density and metallicity of the ambient medium must be assumed in order to compare 
the models with the observations. The metallicity of the gas is an uncertain factor. Most 
spiral galaxies show a decline in metallicity with radius (Ferguson et al. 1998b), (Van Zee 
et al. 1998). Extrapolating the Galactic gradient, with data out to 18 kpc (Smartt & 
Rollestone 1997), to a galactocentric radius of 24 kpc, results in an oxygen abundance 10% 
of the solar value. Similar values were found for outlying H II regions in other galaxies 
(Ferguson et al. 1998b). As oxygen lines contribute significantly to the coohng, the oxygen 
abundance should be a relevant indicator of the total metallicity parameter in the models. 
Therefore, we adopt a metaUicity of 10% of the solar value. 

The scaling relations in Thornton et al. (1998) were used to obtain the model predictions 
listed in Table 2. The model parameters were chosen to cover the relevant range of density 
and metallicity and do not represent a fit to the data. The radius of the model is particularly 
sensitive to the density of the ambient medium. An error of a factor 2 in the density 
translates into a 50% difference in the radius, whereas a factor 10 in metallicity results in a 
20% difference in radius. The mass of the model shells is approximately equally sensitive to 
density and metallicity (20% for a factor 2 variation). The models show a very good overall 
agreement with the observations if the density of the ambient medium is a factor ~ 3 lower 
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than estimated in Section 3.6 {nn — 0.27 cm~^). The age of the model with rin — 0.09 cm~^ 
and log([Z/ZQ]) = —1 is approximately 3 Myr. This is less than the age of GSH 138—01—94 , 
which is consistent with the smaller radius, larger expansion velocity and smaller mass of 
the model. 

The main difference between the models and GSH 138—01—94 is in the mass of the shell, 
which is a factor 5 larger than the model predictions. However, in view of the uncertainty 
in separating emission of GSH 138—01—94 from unrelated Galactic emission, the distance 
to GSH 138—01—94 , the energy release of a supernova explosion adopted in the models, 
the assumed mctallicity, and the ambient density, this is not a major problem. We conclude 
that the observed parameters of GSH 138—01—94 are consistent with those of a very old 
remnant of a single supernova explosion. We note that occasionally supernovae this far from 
the Galactic centre must occur, in view of the evidence for stars as early as spectral type Bl 
in this part of the Galaxy, e.g. De Geus et al. (1993). 

The interpretation of GSH 138—01—94 as a supernova remnant makes it the oldest and 
the largest supernova remnant known. This emphasizes the significance of high-resolution 
H I surveys as a means to find the oldest supernova remnants. It may not be a coincidence 
that such an old but intact supernova remnant is found in the outer galaxy, which is an 
environment very different from that near the solar circle. Individual supernova remnants 
lose their identity when the shock speed drops below the sound speed of the interstellar 
medium, or when they intersect with neighbouring supernova remnants (Ilovaisky et al. 
1974). The timcscalc for these effects arc expected to be longer in the outer Galaxy because 
the density of the interstellar medium is smaller, and supernovae are relatively rare. Also, 
the half- width of the H I disk at 24 kpc is 1.5 kpc, as opposed to 200 pc near the solar circle. 
Therefore, an expanding shell in the outer Galaxy can grow a factor 5 larger before it bursts 
out of the disk and loses its identity as a shell (Mac Low & McGray 1988; Enghsh et al. 
2000). 

Observations of GSH 138—01—94 at wavelengths other than the 21-cm hue of H I would 
be of great value to verify its nature. The significance of the non-detection in the IRAS 
maps is not clear. Emission by dust in the foreground is strong, and far infrared emission 
of the shell may be weak because of a high gas-to-dust ratio and low interstellar radiation 
field in the outer Galaxy. There is no feature at the position of GSH 138—01—94 in the 
34.5 MHz survey of Dwarakanath & Undaya Shankar (1990), or in the 22 MHz survey of 
Roger et al. (1999). However, the upper limit to the flux density of GSH 138—01—94 at 
408 MHz (0.57 Jy) implies a 34.5 MHz flux density less than 7 Jy at 34.5 MHz, assuming a 
spectral index of 1. Therefore, the confusion-hmited sensitivity of these surveys is insufficient 
to expect a detection, and more sensitive observations are needed. X-ray emission may be 
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expected from a hot bubble inside the neutral shell. Figure 7 in Thornton et al. (1998) 
suggests a temperature for the interior that is within a factor 3 of 10^ K. The thermal 
bremsstrahlung spectrum emitted by the hot interior is therefore expected to cut off at a 
photon energy of a few tenths of keV. With a foregroung HI column density 6.3 x 10^^ cm~^, 
these soft X-rays will be heavily obscured. No optical line emission is expected because the 
luminosity of the Thornton et al. (1998) models is less than 0.5% of the peak luminosity 
in this late stage of evolution. Arguably the strongest confirmation of a supernova origin 
of GSH 138—01—94 would be the identification of a central compact object. However, no 
steep-spectrum continuum source could be identified as a pulsar candidate. A source with 
an inverted spectrum was noticed at (Z, 6) = (138?198, — 1?290), but this may well be an 
extragalactic object. Finally, the ongoing arcminute resolution H I surveys of the Galactic 
plane provide the best opportunity to identify similar objects in the far outer Galaxy. 



4.3. Interaction with the environment 

The side of GSH 138—01—94 which faces the Galactic equator is a factor 2 brighter 
than the opposite side. This asymmetry translates directly into a factor 2 in density, consis- 
tent with a density gradient perpendicular to the Galactic equator. This raises the question 
whether the circular shape of GSH 138—01—94 is consistent with this overall asymmetry. 
An estimate of the asymmetry that should be expected can be made by applying the ex- 
pansion law in Thornton et al. (1998), Equation (20). Keeping all parameters constant 
except the ambient density n, we have R ~ Assume that the northern part of 

GSH 138—01—94 expanded in a medium with a density that is on average twice the density 
of the medium on the southern side. Normalizing density and radius to the values in the 
direction of the gradient (south), we have rii = 1, i?i = 1, in the direction of the density 
gradient. Scaling with density, we find n2 — 2, R2 = 0.747 in the opposite direction , and 
ns = 1.5, i?3 = 0.843 for the direction perpendicular to the density gradient. We there- 
fore estimate the axial ratio to be 2i?3/(i?i + R2) = 0.97. The North-South asymmetry in 
brightness is therefore consistent with the nearly circular shape of GSH 138—01—94 . 

The cavities in the upper right quadrant of GSH 138—01—94 discussed in Section 3.4 
may be smaller, presumably younger, shells. Cavity 1 is similar in size to the Cygnus Loop. 
The H II region found by De Geus et al. (1993) is located on the south-eastern edge of Cavity 
1. This is also the region where GSH 138—01—94 appears to make contact with an H I cloud 
(Figure 1). The interpretation of Cavity 1 as a hole in GSH 138—01—94 , and the H II region 
on the perimeter of Cavity 1 is suggestive of progressive star formation. McKee & Ostriker 
(1977) proposed that the passage of a dense cloud through the shell would leave a hole in 
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an expanding shell in the snow plow phase, as if casting a shadow. If this were the case, 
we would expect to observe this cloud. The only candidate, cloud 2 in Digel et al. (1994) 
is much smaller than Cavity 1 (Figure 5). The absence of a radio continuum counterpart of 
both Cavity 1 and Cavity 2 argues against the interpretation that these are young supernova 
remnants. It was argued in Section 4.1 that the non-detection of radio continuum emission 
does not rule out that it is a supernova remnant because of its large age. 

4.4. Implications for galcixies 

We have argued that the outer Galaxy is a less hazardous environment for supernova 
remnants because radiative cooling proceeds at a lower rate, the supernova rate is lower, and 
the scale-height of the gas disk is larger, than in the inner Galaxy. A similar environment 
is usually found in the outskirts of other spiral galaxies, and almost everywhere in dwarf 
irregular and low surface brightness galaxies. GSH 138—01—94 may be the prototype of 
a class of supernova remnant that is representative for such a low metallicity, low density 
environment. The timescale to dissolve completely into the interstellar medium when the 
expansion halts, tdis — Rs/(^v is of the order of 10 Myr. Assuming this process recently 
started in GSH 138—01—94 , at the age of 4.3 Myr, it becomes clear that the shell will 
remain as large density enhancement between 5 and 10 Myr after the supernova explosion. 
We speculate that such old supernova shells may be an important source of dense clouds in 
a low-density environment. This allows star formation to progress for a long time after a 
single trigger event. 

An expanding shell like GSH 138—01—94 would not be resolved in most studies of the 
H I velocity dispersion in galaxies outside the local group, but it would fill a significant part 
of the beam. This may lead to a broader local H I profile, without a counterpart in the 
optical or radio continuum. 

5. Conclusions 

GSH 138—01—94 is a large H I shell at a Galactocentric radius of 24 kpc (heliocentric 
distance 16.6 kpc). Its radius (180 ± 10 pc) and expansion velocity (11.8 ± 0.9 km s~^) 
suggest an age of 4.3 Myr, assuming the expansion law R ~ t^/^. The total mass of the shell 
is 2 X 10^ Mq, and the kinetic energy of the expansion is 3 x 10^^ J. Stellar wind and a 
supernova explosion are considered as possible sources of the shell. 

If GSH 138—01—94 is a stellar wind bubble, the expansion energy of GSH 138—01—94 sets 
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a strong lower limit 60 to the mass of the star if the expansion is driven only by a stellar 
wind. The H II region associated with such a star would have been detected in the radio 
continuum map at 1420 MHz, or optically. There is also no evidence for a young supernova 
remnant inside the H I shell. 

These objections do not exist if GSH 138—01—94 is a supernova remnant. The properties 
of GSH 138—01—94 correspond well with hydrodynamic models of supernova remnants in a 
late stage of evolution. GSH 138—01—94 is then the oldest and largest supernova remnant 
known. Its large size and age may be related to a low-density, low-metallicity environment. 
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Table 2: Properties of GSH 138—01—94 compared with model predictions derived from the 
scahng relations in Thornton et al. (1998). The energy released by the supernova was 
assumed to be 10^^ J. The properties of GSH 138—01—94 are given on the last hne. 





log([Z/Ze]) 


Rs 


1^0 


log(Mtot/M0) 


(cm~^) 




(pc) 


(km s-i) 




0.03 


-1.0 


242.44 


14.36 


4.76 


0.03 


0.0 


192.57 


11.67 


4.49 


0.09 


-1.0 


152.83 


14.20 


4.64 


0.09 


0.0 


121.40 


11.55 


4.37 


0.27 


-1.0 


96.34 


14.05 


4.53 


0.27 


0.0 


76.53 


11.42 


4.26 


0.27 


-l.O'^ 


180±10 


11.8±0.9 


5.3±0.3 



'^Metallicity estimated from the Galactic oxygen abundance gradient (Smartt & Rollestone 1997). 
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